AD-A159  682  INVESTIGATION  OF  THE  SPATIALLV  ISOTROPIC  CONPONENT  OF  1/1 

THE  LATERALLV  AVERA.  .  <U)  CHICAGO  UNIV  IL  JANES  FRANCK 
INST  C  F  VU  ET  AL.  AUG  85  TR-6  N88814-77-C-B248 

F/G  28/18 


UNCLASSIFIED 


NL 


•  •  ] 


unclassifi ed 

•  > 


SICUWTY  CLASSIFICATION  of  this  PAGE  (Whmt  Dmtm  Entmd) 


RETORT  DOCUMENTATION  RAGE _  aEpSt^SSXISg'VoRM 


2.  GOVT  ACCESSION  NOJ  1.  RECIPIENT'S  CATALOG  NUMEER 


'  INV^'riGAff^  OF  TliE  SPATIALLY  ISOTROPIC  COMPONENT  *'  TYPe  °r  RePO,,T  *  ”"'°0  covered 
OF  THE  LATERALLY  AVERAGED  MOLECULAR  IIYDROGEN/Ag(  111.)  Interim  Technical  Report 

PlIYS  I  SORPTION  POTENTIAL  _ _ _ _ 

•  performing  oro.  report  number 


I.  CONTRACT  or  grant  NUMBERS 

N00014-77-C-0240 


autmort*; 

Chien-fan  Yu,  K.  Birgitta  Whaley,  C.  S.  Hogg  and 
Steven  J.  Sibener 


performing  organization  name  ano  aooress 

The  University  of  Chicago 

The  James  Franck  Institute 

5640  South  Ellis  Avenue,  Chicago,  IL  60637 


II.  CONTROLLING  OFFICE  NAME  ANO  AOORESS 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  0  WORK  UNIT  NUMBERS 


NR- 392-023 


IS.  REPORT  OATS 


The  Office  of  Naval  Research 
Physical  Sciences  Division  (Code  421) 

800  North  Quincy  Str. ,  Arlington,  VA  22217 


’  NAME  A  AODRESSfJf  dj/termi  tnm  Controlling  Ol  tic,)  IS.  SECURITY  CLASS.  (•!  (Ala 

unclassified 


August  1985 


IS.  NUMBER  OF  PAGES 

70 


ISa.  DECLASSIFICATION/  DOWNGRADING 
SCHEDULE 


IG.  DISTRIBUTION  STATEMENT  (ol  till,  Roport) 


This  document  has  been  approved  for  public  release  and  sale;  its 
distribution  is  unlimited. 


DTIC 


IB.  SUPPLEMENTARY  notes 


1 

Prepared  for  Publication  in  Journal  of  Chemical  Physics 

lr  '  cl 


WORDS  (CmiIum  on  roooroo  otfdo  //  fioooooorr  a ntf  Identity  by  block  »uobwj 

Gas-surface  interaction  potentials;  physisorption^>s elective  adsorption; 
para-hydrogen  generation;  gas-surface  scattering  resonances^^ 


ABSTRACT  fCanllDM  an  raaaraa  <l4a  If  nacaaaarr  aN  Manflfr  by  UmI  nombM) 

A  comprehensive  study  of  the  spatially  isotropic  component  of  the  laterally 
averaged  molecular  hydrogen/Ag( 1 1 1)  physisorption  potential  is  presented.  Dif¬ 
fractive  selective  adsorption  scattering  resonances  for  rotationally  state- 
selected  1  i 2  and  D2on  Ag(lll)  have  been  mapped  out  as  a  function  of  incident  polar 
angle  for  several  crystal  azimuths  and  beam  energies.  These  resonances  have  been 
used  to  determine  the  bound  eigenvalues,  and  subsequently  the  shape,  of  the  poten¬ 
tial  well.  Best  fit  Lennard-Jones ,  Morse,  variable  exponent,  and  exponential-3 
potentials  having  well  depths  ^32  meV  are  derived  from  the  data.  These  measure- 


00  ,  '2Tn  1473 


EDITION  OF  I  NOV  ••  IS  OBSOLETE 
S/N  0.02-  LP- 014*4401  > 


unclassified 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  fWMn  Oat*  BHIarW) 


*5  10  02  020 


unclass ified 

tCCUM TV  CLASSIFICATION 


II! 


ments  are  supported  by  rotationally  inelastic  scattering  measurements  for  HD  and 
exact  close-coupled  quantum  scattering  calculations.  Debye-Waller  attenuation 
measurements  are  also  presented  for  H2,  D2,  and  HD.  The  ability  to  detect  these 
diffractively  coupled  resonances  on  a  closest-packed  metallic  surface,  i.e.  a 
surface  of  extremely  low  corrugation,  suggests  that  such  measurements  can  be 
carried  out  on  a  much  wider  class  of  surfaces  than  previously  envisioned. 


A  COMPLETE  COPY  OF  THE  MANUSCRIPT  WILL  BE  SENT  UPON  REQUEST. 


unclassified 

•■CUN! TV  CLASSIFICATION  OF  THIS  FAOtfWfcaa  Data  SMara4 

I  i 


S/N  0103- l*  014-4401 


INVESTIGATION  OF  THE  SPATIALLY  ISOTROPIC  COMPONENT 


OF  THE  LATERALLY  AVERAGED  MOLECULAR 
HYDROGEN/Ag( 111)  PHYSISORPTION  POTENTIAL 

Chien-fan  Yi/a^,  K.  Birgitta  Whaley^),  C.  S.  Hogg 
and  Steven  J.  Sibener^ 

Department  of  Chemistry  and 
The  James  Franck  Institute 
The  University  of  Chicago 
5640  South  Ellis  Avenue 
Chicago,  Illinois  60637 


In  Press,  Journal  of  Chemical  Physics 

^Alfred  P.  Sloan  Fellow 
PACS:  79.20 .Rf 

(a)  Current  Address:  Dept,  of  Chemistry,  Columbia  University,  New  York 

(b)  Current  Address:  Dept,  of  Chemistry,  Tel  Aviv  University,  Israel 


ABSTRACT 


l 

V 

A  comprehensive  study  of  the  spatially  isotropic  component  of  the  laterally 
averaged  molecular  hydrogen/Ag( 111)  physisorption  potential  is  presented.  Dif¬ 
fractive  selective  adsorption  scattering  resonances  for  rotationally  state-selected 
and  on  Ag(lll)  have  been  mapped  out  as  a  function  of  incident  polar  angle 
for  several  crystal  azimuths  and  beam  energies.  These  resonances  have  been  used 
to  determine  the  bound  eigenvalues,  and  subsequently  the  shape,  of  the  potential 
well.  Best  fit  Lennard- Jones ,  M&csje,  variable  exponent,  and  exponential-3  poten- 
tials  having  well  depths  of~^32  meV  are  derived  from  the  data.  These  measurements 
are  supported  by  rotationally  inelastic  scattering  measurements  for  HD  and  exact 
close-coupled  quantum  scattering  calculations.  Debye-Waller  attenuation  measure¬ 
ments  are  also  presented  for  H^,  D^,  and  HD.  The  ability  to  detect  these  diffrac- 
tively  coupled  resonances  on  a  closest-packed  metallic  surface,  i.e.  a  surface 
of  extremely  low  corrugation,  suggests  that  such  measurements  can  be  carried  out 
on  a  much  wider  class  of  surfaces  than  previously  envisioned.  Qr’\  c^'s  V'  u  l 
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Introduction 

There  is  currently  great  interest  in  formulating  an  accurate  description  of 
gas-surface  interaction  potentials  for  a  wide  variety  of  systems.  This  is  clearly 
indicated  by  the  large  body  of  recent  experimental  and  theoretical  work  in  this 
area  [1,  2].  The  motivation  for  these  studies  is  rather  straightforward:  know¬ 
ledge  of  the  interaction  potential  is  a  prerequisite  for  a  complete  understanding 
of  many  important  heterogeneous  processes  such  as  oxidation,  catalysis,  condensa¬ 
tion,  and  energy  transfer,  and  is  intimately  related  to  the  structure,  charge 
density  distribution,  and  dynamics  of  the  surface.  The  central  theme  of  this 
paper  will  be  to  present  a  detailed  study  of  the  spatially  isotropic  component 
of  the  laterally  averaged  molecular  hydrogen/ Ag( 111)  physisorption  potential. 

This  will  be  accomplished,  in  part,  by  using  diffractive  selective  adsorption 
(DSA)  scattering  resonances  to  determine  the  bound  eigenvalues  of  the  potential. 
Resonances  of  this  type  were  first  observed  by  Estermann  and  Stern  in  1930  [3], 
and  subsequently  explained  by  Lennard-Jones  and  Devonshire  [4].  This  phenomenon, 
due  to  both  the  excellent  energy  resolution  of  contemporary  gas-surface  scatter¬ 
ing  instruments  and  the  implementation  of  accurate  quantum  scattering  calculations, 
provides  the  most  valuable  means  of  determining  laterally  averaged  physisorption 
potentials  at  this  time. 

For  the  physisorption  interaction  of  atoms  and  molecules  with  surfaces  of 
perfect  periodicity,  the  interaction  potential  with  a  stationary  lattice  can  be 
expanded  as  a  Fourier  sum  over  the  surface  reciprocal  lattice  vectors,  6.  For 
atom-surface  systems 

V(.z,  R)  =  I  V_(z)  e1®’?  (1.1) 

G  “ 

with  z  the  perpendicular  distance  of  the  atom  from  the  surface  and  R  the  surface 


by  Vg(z,  6,  4>,  r)  where  9,  <j>,  and  r  are  the  additional  spherical  polar  internal 
coordinates  of  the  diatom  and  z,  R  now  refer  to  the  location  of  the  molecule 
center  of  mass  in  the  surface  fixed  coordinate  system.  Neglecting  in  the  first 
approximation  the  internal  degrees  of  freedom  for  the  molecular  case,  the  later¬ 
ally  averaged  interaction  potential  VQg(z)  can  be  divided  into  two  components: 
an  attractive  long  range  polarization  potential,  Vgg  ^(z)  and  a  short  range  re¬ 
pulsive  potential,  VQQ  R(z).  The  asymptotic  form  of  the  van  der  Waals  term  is 
well  established  as  [5,  6] 

Wz> '  -c3/(z-y3  (i-2> 

where  z  is  the  normal  distance  between  the  adparticle  center  of  mass  and  the 
surface.  Both  and  the  position  of  the  reference  plane,  Zq,  have  been  calculated 
for  a  number  of  atom-solid  systems  [6].  Ab  initio  calculation  of  the  repulsive 
component  requires  knowledge  of  both  the  electronic  structure  of  the  surface  and 
the  adparticle,  and  inclusion  of  both  Coulombic  and  exchange  terms.  Several  cal¬ 
culations  have  been  made  for  the  interaction  of  helium  with  jellium  surfaces  [7]. 
Recently  Harris  and  Liebsch  [8]  have  calculated  a  repulsive  potential  for  He/Cu(110) 
containing  two  adjustable  parameters  which  gives  good  agreement  with  the  available 
scattering  data  on  this  system  [9].  Since  their  approach  applies  to  any  adparticle 
with  a  closed  shell  electronic  structure,  it  was  also  possible  to  provide  an 
estimate  of  the  corresponding  H2/Cu(110)  potential  [10],  including  an  estimate  of 
the  anisotropic  contribution  to  the  laterally  averaged  physisorption  potential 
Vqo(z»  6.  4> »  r) .  Such  ab  initio  calculations  on  real  systems  present  considerable 
difficulties.  In  this  respect,  a  recent  suggestion  that  the  repulsive  He-surface 
interaction  potential  can  be  approximated  as  a  local  function  of  the  surface 
electronic  charge  density  [11]  may  indicate  a  viable  practical  alternative 
approach  to  the  ab  initio  calculation  of  other  physisorption  potentials  [12]. 

By  comparison,  the  approach  we  take  here  is  the  empirical  determination  of  the 
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laterally  averaged  potential  which  is  the  result  of  extensive  molecule- 
surface  scattering  experiments  and  theoretical  scattering  calculations. 

Molecular  and  atomic  beam  scattering  experiments  allow  determination  of 
the  physisorption  well  depth  D,  the  hound  states  of  the  laterally  averaged 
potential  (and,  therefore,  the  shape  of  the  well),  but  not  its  absolute 
position  with  respect  to  the  surface  atoms.  In  addition,  molecular  scat¬ 
tering  gives  information  on  the  anisotropic  (i.e.  orientation  dependent) 
component  of  the  laterally  averaged  potential.  Extensive  experiments  and 
calculations  on  scattering  of  the  hydrogen  isotopes  Hg>  Dg,  and  HD  from 
Ag(lll)  have  allowed  us  to  make  a  detailed  analysis  of  the  full  laterally 
averaged  physisorption  potential  on  this  closely  packed  and  extremely  smooth 
surface.  In  this  paper  we  present  the  results  relating  to  the  spatially 
isotropic  component  of  the  molecular  hydrogen/Ag( 111)  potential.  The  aniso¬ 
tropic  component  will  be  presented  in  a  forthcoming  paper  [13]. 

In  these  experiments  we  measure  both  diffractive  selective  adsorption,  DSA, 
and  rotationally  mediated  selective  adsorption,  RMSA,  resonances  for  and  on 
Ag(lll).  This  is  one  of  the  first  studies  reported  in  which  DSA  techniques  have 
been  successfully  used  on  a  closest  packed  (i.e,,  extremely  low  corrugation)  metallic 
surface  to  extract  information  on  a  physisorption  potential,  and  demonstrates,  to¬ 
gether  with  data  on  the  hydroqen/Pt(l 11)  system  [14],  that  such  measurements  can  be 
carried  out  on  a  much  wider  class  of  surfaces  than  previously  thought  possible. 
Experiments  in  this  field  have  focussed  primarily  on  studies  of  DSA  resonances  in 
the  scattering  of  light  particles  (He,  H  and  H 2)  from  highly  corrugated  surfaces 
of  dielectrics  like  Li F  [15],  NaF  [16],  graphite  [17-19],  and  stepped  metal  surfaces 
of  copper  [9].  The  scarcity  of  DSA  results  on  metals  may  be  attributed  to  the  very 
low  non-zero  order  diffraction  probabilities  resulting  from  the  small  corrugation 
of  metal  surfaces  [20].  However,  the  observation  that  H^  diffraction  peaks 
on  Ag(lll)  are  an  order  of  magnitude  stronger  than  He  peaks  [21]  under  equiva¬ 
lent  incident  conditions  suggests  that  scattering  from  Ag( 111)  is  a  promising 
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candidate  for  molecular  DSA  studies.  This  observation  was  explained  by  Hill 
et  al .  [22]  in  terms  of  a  lateral  variation  of  the  van  der  Waals  interaction, 
and  more  recently  by  Liebsch  and  Harris  [23]  in  terms  of  differences  in  the 
repulsive  component  of  the  interaction.  Silver  is  a  good  substrate  for  DSA 
experiments  since  it  is  chemically  inert  with  respect  to  hydrogen  chemisorption, 
so  the  beam  scattering  experiments  can  be  done  at  surface  temperatures  down  to 
^30  K  without  interference  due  to  hydrogen  loading  of  the  surface.  The  use  of 
low  surface  temperatures  together  with  molecular  beams  of  very  low  energy  mini¬ 
mizes  Debye-Waller  type  attenuation,  and  enables  us  to  overcome  the  difficulty 
of  detecting  resonances  coupled  via  very  low  diffraction  probabilities.  Low  cor¬ 
rugation  actually  facilitates  the  analysis  of  DSA  experiments  since  it  virtually 
eliminates  any  complications  due  to  band  structure  [24]. 

Past  DSA  studies  on  molecular  hydrogen  have  treated  it  as  a  structureless 
particle  [18].  However,  it  is  to  be  expected  that  bound  state  energy  levels 
obtained  from  selective  adsorption  measurements  will  be  affected  by  the  spatially 
anisotropic  component  of  the  molecule-surface  interaction  potential.  Recent  low 
energy  electron  energy  loss  (EELS)  experiments  on  the  physisorption  of  H^  and  D 2 
on  Ag(lll)  [25]  and  of  H2>  D2  and  HD  on  Cu(110)  [26]  indicate  that  molecular 
hydrogen  behaves  as  a  nearly  free  3-dimensional  rotor  on  these  surfaces.  However, 
the  resolution  was  too  low  (12  meV  for  Ag(lll),  3  meV  for  Cu(100))  for  quantitative 
measurement  of  any  deviation  from  the  unhindered  results.  In  EELS  experiments 
the  spatial  anisotropy  should  manifest  itself  in  slight  shifts  of  the  observed 
vibration-rotation  transitions  for  physisorbed  |J,m>  states,  and  as  a  consequence 
of  this,  in  linewidth  changes  for  spectra  obtained  as  a  convolution  over  various 
m  sub-states.  Calculations  on  the  RMSA  resonances  of  HD  on  Pt ( 111)  and  HD  on 
Ag (111)  have  shown  that  RMSA  resonances  of  physisorbed  HD  are  consistent  with 
bound  states  of  the  isotropic  component  plus  a  first  order  level  shift  deriving 
from  the  anisotropic  component  of  the  interaction  potential  [27,  28]. 


The  inadequacy  of  a  free  rotor  approximation  is  also  evidenced  in  other  calcu¬ 
lations  for  the  HD/Pt( 111)  system  [29].  This  shift  in  general  will  be  J 
and  m  dependent.  Similar  |J,m>  dependent  shifts  should  also  appear  in  the  bound 
state  spectrum  of  the  physisorption  well  when  probed  by  selective  adsorption 
scattering  resonances.  Comparison  of  OSA  resonance  energies  for  n-H^  (J  =  1: 

J  =  0  is  3:1)  and  p-H^  and  both  DSA  and  RMSA  resonances  for  n-D^  (J  =  1:  J  =  0 
is  1:2)  and  0-D2  have  enabled  us  to  directly  detect  an  m  state  averaged  shift 
in  these  experiments. 

The  plan  for  the  remainder  of  the  paper  is  as  follows:  Sections  2,  3  and 
4  describe  the  apparatus,  construction  of  the  p-H2  and  o-D^  converter,  and  the 
experimental  procedure.  Section  5  describes  the  Debye-Waller  analysis  and  pro¬ 
blems  arising  from  this  for  the  molecular  case.  The  kinematics  of  selective 
adsorption  resonances  are  analyzed  in  detail  in  Section  6  and  the  connection 
between  DSA  and  RMSA  is  given  explicitly.  The  experimental  selective  adsorption 
data  are  presented  in  Section  7  and  quantum  number  assignments  discussed  in 
Section  8.  The  isotopic  study  of  p-H2  and  o-D2  (J  =  0)  resonances  according  to 
the  semiclassical  method  of  Le  Roy  [30]  gives  two  possible  quantum  number  assign¬ 
ments:  we  distinguished  between  these  by  performing  rotationally  inelastic 
quantum  scattering  calculations  for  HD  on  Ag (111)  for  model  potentials  derived 
from  each  assignment,  and  comparing  the  results  to  experimental  HD  rotationally 
inelastic  transition  probabilities.  Section  9  contains  a  brief  comparison  with 
a  recently  proposed  scaling  relation  of  Vi  dal i  et  al .  [31],  which  further  dis¬ 
tinguishes  the  two  quantum  assignments.  Thermal  desorption  experiments  are 
described  in  Section  10.  Finally,  Section  11  presents  a  brief  summary  and  our 
conclusions . 
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II .  Apparatus 

The  scattering  apparatus  used  in  these  experiments  has  been  previously 

described  [28,  32].  It  consists  of  four  major  sections:  the  beam  source  and 

vacuum  manifold,  rotating  detector,  crystal  assembly  and  associated  surface 

analysis  instrumentation,  and  data  acquisition  electronics. 

The  nozzle  beam  is  differentially  pumped  in  three  stages  in  a  non-bakeable 

-3  -5 

source  chamber,  with  typical  operating  pressure  being  5  x  10  ,  1  x  10  ,  and 

5  x  Kf7  torr  for  the  regions  before  the  skimmer,  before  the  beam  collimator, 

and  containing  the  chopper,  respectively.  The  main  UHV  scattering  chamber  is 

separated  from  the  three  source  chamber  regions  by  a  bakeable  buffer  region 

_g 

whose  pressure  rises  to  5  x  10  torr  when  the  beam  is  on.  This  buffer  section 

is  separated  from  the  chopper  section  by  a  sliding  gate  valve  which  contains  the 

0.5  mm  diam.  beam  defining  aperature.  Supported  by  a  400  1/s  ion  pump  and  a 

titanium  sublimation  pump,  the  crystal  chamber  has  a  base  pressure  of  7  x  10”11 

-9 

torr,  which  rises  as  high  as  1  x  10  torr  when  the  beam  is  on.  The  incident 
beam  is  collimated  to  an  angular  divergence  of  0.1°,  which  corresponds  to  a  0.5mm 
diam.  beam  spot  at  the  crystal. 

Particle  detection  is  accomplished  with  a  doubly  differentially  pumped  quad¬ 
ruple  mass  spectrometer,  which  rotates  about  a  horizontal  axis  in  the  plane  con¬ 
taining  the  incident  beam  and  the  surface  normal.  The  electron  bombardment  ionizer 
is  located  14.45  cm  from  the  sample,  and  is  collimated  to  view  the  crystal  with 
0.67°  resolution.  The  detector  can  be  rotated  manually  or  under  computer  control. 

The  crystal  manipulator  provides  x,  y,  and  z  translational  degrees  of  freedom, 
incident  polar  angle  and  surface  azimuthal  angle  scanning,  and  in-plane  and  out-of¬ 
plane  tilt  adjustment.  The  retractable  CMA  Auger  spectrometer  with  a  side  mounted 
ion-sputtering  gun  can  be  placed  under  or  away  from  the  sample  before  and  after 
ion  sputtering  or  taking  Auger  spectra. 
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Data  acquisition  and  instrument  control  is  carried  out  with  a  PDP-11  mini¬ 
computer  which  is  interfaced  to  the  scattering  apparatus  via  CAMAC.  Time-of- 
flight  measurements  are  made  with  a  255  channel  scaler  which  has  a  minimum  dwell 
time  of  0.25ysec/channel .  Molecular  beam  gating  is  carried  out  with  a  15.24  cm 
diam.  chopper  which  contains  two  narrow  slots  for  velocity  analysis  (0.5mm  wide) 
and  two  50%  duty  cycle  patterns  for  angular  distribution  measurements.  Angular 
distributions  are  obtained  with  two  32  bit  scalers  which  are  synchronized  to  the 
chopper  in  order  to  collect  signal  plus  background  (slot  open)  and  background 
only  (slot  closed)  information.  DSA  signals  were  collected  by  a  Keithley  elec¬ 
trometer,  demodulated  by  a  lock-in  amplifier,  passed  to  a  digital  voltmeter,  and 
finally  read  into  the  computer. 

III.  p-H2  and  o-D2  Generation 

The  molecular  nuclear  state  degeneracies  of  H 2  and  D?  give  rise  to  the  para/ 
ortho  population  ratio  of  1:3  for  (normal-H2)  and  1:2  for  D2  (normal-D2)  in 
the  high  temperature  limit.  These  ratios  will  not  change  in  the  course  of  super¬ 
sonic  jet  expansion.  For  an  ideal  expansion  the  terminal  rotational  state  distri¬ 
bution  of  H2  will  be  25%  J  =  0  state  (p-H2)  and  75%  J  =  1  state  (o-H2).  Similarly 
the  D2  beam  will  be  33%  J  =  1  ( p-D^ )  and  about  67%  J  =  0  (o-D2).  In  order  to  pre¬ 
pare  pure  rotational  ground  state  beams  of  H2  and  D2,  conversion  to  p-H2  and  o-D2 
has  to  be  made  before  the  jet  expansion. 

An  in-line  p-hydrogen/o-deuteri urn  converter  was  constructed  as  shown  in 
Figure  1.  The  converter  consists  of  a  coaxial  stainless  steel  tube  arrangement 
which  contains  a  nickel  salt-based  catalyst,  APPACHI  [33],  in  its  outer  tube.  A 
sintered  stainless  steel  filter  element  was  welded  to  the  end  of  the  inner  tube 
to  prevent  the  catalyst  from  entering  this  region  of  the  converter.  The  catalyst 
was  activated  by  heating  at  150  to  175°  C  for  a  few  hours  while  flowing  a  small 
stream  of  hydrogen  gas  through  it  [34],  In  operation,  the  converter  was  immersed 


in  a  liquid  hydrogen  dewar  to  maintain  a  very  low  temperature,  ensuring  excellent 
conversion  to  J  =  0  hydrogen. 

When  n-H^  or  n-D^  are  brought  into  the  converter,  exposure  to  the  strong 
magnetic  field  gradients  at  the  surface  of  the  catalyst  will  quickly  promote 
equilibration  at  this  low  temperature,  i.e.  99.8%  p-H^  and  97.8  o-D^.  The  almost 
pure  p-H^  or  o-D^  then  exits  through  the  inner  tube  to  the  nozzle.  The  converter 
can  operature  at  pressures  up  to  2000  psig,  more  than  that  required  for  producing 
for  a  good  jet  expansion.  A  by-pass  line  was  also  incorporated  to  facilitate 
immediate  switching  between  n-/p-H2  or  n-y/o^. 

Comparison  of  the  rotationally  inelastic  scattering  transition  probabilities 
for  H2  or  D2  with  and  without  conversion  can  be  used  to  check  the  effectiveness  of 
the  converter.  Figures  2  and  3  show  the  angular  distributions  of  and  scat¬ 
tering  from  Ag(lll)  with  or  without  conversion.  By  subtracting  the  diffuse  back¬ 
ground  and  deconvoluting  the  instrumental  factors  from  the  specular  and  inelastic 
peaks,  relative  probabilities  of  the  J  =  0  -*  2  transition  to  that  of  the  specular 
beam  can  be  obtained,  and  are  listed  in  the  last  column  of  Table  1. 

The  rotational  state  populations  in  the  expanded  jet  can  be  estimated  by  solv¬ 
ing  the  energy  balance  equation  [35]: 


!  Vo  *  MV  --W  Vp  *  MV 


(3.1) 


where  Tq  is  the  source  temperature,  u  the  stream  velocity,  Tp  parallel  temperature 
of  the  beam,  and  T^  the  rotational  temperature  of  the  beam.  A  Boltzmann  distri¬ 
bution  of  rotational  states  is  assumed  in  the  non-commuting  isothermal  canonical 
ensembles  of  para  and  ortho  components.  The  results  are  also  shown  in  Table  1. 

It  is  clear  that  the  population  of  the  J  =  0  ^  increases  by  a  factor  of  four 
after  conversion  which  is  close  to  the  ratio  of  the  relative  probability  of  the 
J  =  0  -*■  2  transition  of  p-f^  to  that  of  n-h^.  Similarly,  the  J  =  0  D2  increases 
by  a  factor  of  1.5  which  is  approximately  the  ratio  of  the  relative  probability  of 


T- 
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fall  on  three  sets  of  concentric  arcs  centered  around  G  =  (10),  (01),  and  (II). 
state  energies  for  J  =  0  can  be  directly  calculated  from  arcs  using  Eq.  (6.7) 
to  analyze  the  data.  Similar  measurements  have  also  been  made  using  o^,  as 
shown  in  Figures  13  and  14.  In  addition  to  the  0-D2  DSA  dips,  two  RMSA  dips 
also  appear  (dash-dot-dot  lines)  as  broader,  azimuth-independent  features,  for 
which  bound  state  energies  en  jQ  can  be  calculated  using  Eq.  (6.15).  Another 
observation  is  that  no  discernible  splittings  occur  at  <)>  =  30°  between  the 
resonances  coupled  by  the  (11)  and  (10)  lattice  points,  i .e.  the  two  groups  of 
arcs  coalesce  smoothly  at  this  azimuth.  This  absence  of  band  structure  at  the 
crossing  points  of  the  arcs  is  presumably  due  to  the  very  small  non-zero  order 
Fourier  components  of  the  periodic  surface  potential  [24]. 

Figures  15  and  16  show  the  specular  intensity  spectra  for  n-H2  and  n-D^, 
respectively,  taken  at  three  different  beam  energies  along  the  <112>  direction. 
The  resonances  can  be  seen  to  get  stronger  and  broader  as  the  beam  energy  is 
lowered.  These  energy  dependent  runs  serve  to  confirm  our  level  assignments, 
and  demonstrate  that  low  energy  DSA  spectra  are  more  sensitive,  and  have  better 
energy  resolution,  than  higher  energy  experiments.  The  low  energy  data  has  in 
fact  enabled  detection  of  an  extra  bound  state  near  the  top  of  the  well  for  both 
and  D^ .  The  estimated  accuracy  of  the  bound  state  resonance  energies  deter¬ 
mined  with  the  low  energy  incident  beam  are  ca.  0.2  meV  for  and  0.15  meV  for 

The  selective  adsorption  loci  for  both  isotopes  have  also  been  mapped  out 

in  the  K  -K  reciprocal  lattice  plane  at  beam  energies  of  approximately  35  and 
x  y 

80  meV,  yielding  results  in  agreement  with  the  low  energy  data.  The  higher 
energy  spectra  contain  resonance  dips  due  to  an  additional  bound  state, 
cQ  %  -25.7  meV. 


Bound 
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were  achieved  by  redesigning  the  cryogenic  connections  to  the  crystal  mount, 
thereby  decreasing  Oebye-Waller  attenuation  of  the  specular  beam.  The  crystal 
mount  was  also  modified  to  include  a  pointer  and  vernier  scale  which  allowed 
the  azimuthal  angle  to  be  accurately  set  for  orientations  falling  between  the 
<ll2>  and  <10l>  symmetry  directions  (see  Section  IV).  In  addition,  a  new  beam 
source  was  constructed  in  which  the  nozzle  assembly  was  directly  mounted  on  a 
liquid  cryogen  reservoir  located  in  the  nozzle-skimmer  vacuum  region.  This 
allowed  us  to  produce  an  extremely  stable  H2  beam  of  lower  energy  than  was  pre¬ 
viously  available  in  our  laboratory.  The  actual  beam  pumping  manifold  was  also 
rebuilt,  eliminating  booster  pumps  in  the  nozzle-skimmer  region,  and  replacing 
them  with  high  throughput  diffusion  pumps.  This  last  modification  resulted  in 
the  production  of  more  intense  hydrogen  beams  of  narrower  velocity  dispersion. 
Finally,  in  order  to  address  concerns  pertaining  to  the  then  unknown  magnitude 
of  the  resonance  energy  shifts  occurring  due  to  the  J  f  0  components  in  the  in¬ 
cident  beam,  an  in-line  p-H2/o-D2  converter  was  designed  and  fabricated  (Section 
III). 

Figure  11  displays  eight  p-H^  resonance  spectra  taken  at  successive  azimu¬ 
thal  angles.  DSA  dips  appearing  in  these  runs  are  linked  by  distinguishable  lines 
which  group  the  resonances  according  to  the  pertinent  coupling  reciprocal  lattice 
vector  (mn).  Individual  dips  are  designated  by  (mn\,  where  vis  the  quantum  num¬ 
ber  of  the  level.  (The  manner  in  which  these  assignments  are  made  will  be  dis¬ 
cussed  in  detail  in  Section  VIII.)  The  positions  of  the  dips  attributable  to  a 
given  coupling  G  vector  clearly  shift  in  a  similar  manner  as  the  crystal  azimuth 
is  varied.  This  feature  is  readily  seen  in  Figure  12,  which  is  the  corresponding 
plot  for  the  DSA  resonances  presented  in  the  surface  reciprocal  lattice  Kx'Ky 
plane.  The  resonance  dip  positions,  represented  by  points  in  this  figure,  clearly 
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Voo(z’  6)  =  vSo(z)  +  voo(z)  P2(cose)  (6*21) 

and 

£n!jm  =  <n I voo^ * n><Jm I p2^ cos0) I Jm>  (6.22) 

2 

Even  though  the  strength  of  the  anisotropic  term  vQo(z)  is  very  weak  for  H 2 

and  D2>  the  effects  of  these  shifts  are  large  enough  to  be  observed  by  comparison 

of  the  resonance  positions  and  widths  for  normal  H2  (J  =  1  :  J  =  0  is  3:1)  and 

normal  D2  (J  =  1  :  J  =0  is  1:2)  with  those  of  para-H2  and  ortho-D2  respectively. 

A  preliminary  report  of  these  shifts  has  already  appeared  [51].  In  addition, 

since  for  J  =  J.  =  0  (m  =  m.  =0)  is  zero,  comparison  of  the  DSA  resonance 

i  i  n , Jm 

energy  with  the  RMSA  resonance  energy  from  the  same  level  n  when  J..  =  0  provides 

2 

a  sensitive  means  of  isolating  the  effect  of  the  anisotropic  component  vqo(z), 
provided  the  isotropic  component  v°Q(z)  is  known.  This  will  be  discussed  in 
detail  in  a  forthcoming  paper  [13].  We  now  present  the  selective  adsorption 
results  for  normal  and  state  selected  (J  =  0)  H2  and  D2- 

VII .  Selective  Adsorption  Resonances:  Experimental  Results 

During  the  early  stages  of  this  work  only  n-H2  and  n-D2  were  available  for 
use  in  searching  for  the  presence  of  DSA  resonances  on  Ag (111).  The  existence  of 
discernible  DSA  dips  was  first  confirmed  by  comparing  specular  intensity  data, 
taken  as  a  function  of  incident  polar  angle,  at  azimuthal  angles  symmetric  about 
the  <ll2>  direction,  as  shown  in  Figure  10  for  4>  =  ±  20°.  In  this  figure  weak  but 
distinctive  pairs  of  dips  can  be  seen  in  the  4>  =  20°  and  <}>  =  -20°  spectra,  as  in¬ 
dicated  by  the  six  dashed  lines.  The  presence  of  these  low  sianal-to-noise  fea¬ 
tures  in  both  spectra,  as  required  by  the  symmetry  of  the  reciprocal  lattice  with 
respect  to  the  <112>  direction,  provided  us  with  the  necessary  motivation  to  map 

out  these  resonances  across  the  K  -K  plane. 

x  y 

Several  instrumental  improvements  were  then  made  with  the  intent  of 
increasing  the  quality  and  precision  of  the  data.  Colder  surface  temperatures 


For  the  rotationally  elastic  case,  J  =  J^,  this  reduces  to  pure  DSA,  with 


resonance  condition 

s2i 


ki  =  JO)  +  J1)  (K  +  g)2 

£n  en,Jm  2M  '5i 


(6.20) 


and  we  see  that  the  DSA  resonance  energies  also  show  first  order  shifts 

An  important  distinction  between  pure  DSA  and  pure  RMSA  resonances  is  that 
the  former  are  azimuthal ly  dependent,  Eq.  (6. 20 X  while  the  later  show  no  azimu¬ 
thal  dependence,  Eqs .  (6.14-15).  Mixed  DSA-RMSA  resonances  will  also  be  azimuthally 
dependent,  Eqs .  (6.18-19).  No  evidence  for  such  azimuthally  dependent  resonances 
have  been  observed  to  date  in  HD  (specular)  rotationally  inelastic  scattering 
from  surfaces.  The  relative  magnitude  of  the  rotational  coupling  on  smooth  metal 
surfaces  is  so  much  greater  than  diffractive  coupling  for  HD,  that  the  small 
azimuthally  dependent  variations  in  the  specular  beam  rotational  distributions 
due  to  diffractive  coupling  are  not  detectable  at  the  present  sensitivity. 

However,  the  situation  is  very  different  for  the  homonuclear  hydrogen 
isotopes.  Although  again  no  mixed  DSA-RMSA  resonances  have  been  observed  yet, 
the  very  much  smaller  rotational  inelasticity  makes  rotational  and  diffraction 
excitation  probabilities  of  comparable  magnitude  and  we  find  resonances  in  the 
specular  beam  which  fall  into  two  groups.  These  are  i)  azimuthally  dependent 
resonances,  due  to  DSA  (Eq.  6.20)  and  ii)  azimuthally  independent  resonances, 
due  to  pure  RMSA  (Eqs.  (6.14,  6.15)).  The  first  order  level  shift  in  these 
resonance  conditions  thus  implies  a  dependence  on  the  J  and  m  quantum  numbers  of 
the  initial  state  for  DSA  resonance  energies,  and  a  dependence  on  J  and  m  of 
both  the  initial  state  and  the  closed  rotational  channel  coupled  to  this  at  the 
resonance  for  RMSA.  For  and  the  laterally  averaged  potential  may  be 
represented  by  the  first  two  symmetric  terms  in  a  Legendre  expansion 
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where  v“0(z)  is  the  spatially  isotropic  component  which  is  the  same  for  all  iso¬ 
topes  ( H2 »  D 2 >  and  HD).  The  first  order  approximation  to  the  rotational  resonance 
energies  is  then  given  by  eigenvalues  of  the  rotationally  decoupled  equations 


$  +  kz  +  T  VJi  +  1)  -  T  J(J  ♦  »)  -  J  »0°0<Z>1W2»  ■ 


*  =  “<2)  <Jmlpa‘cos0>lJm>Wz> 


(6.13) 


For  the  homonuclear  isotopes,  vQ0(z)  is  equal  to  zero  [13].  Provided  v°Q(z)  is 

~  0 

greater  than  vqo(z),  the  first  order  approximation  to  the  eigenvalues  en  jm  of 
Eq.  (6.11)  are  then  given  by 


=  ♦  d 

n,Jm  z  21  i'  i  ' 


(6.14) 


en  +  21  J(J  +  1)  +  Z  ^mlP^Ccose)  |Jmxn  jv^Q(z)  |n>  (6.15a) 


-  e(°)  +  ed) 

n,Jm  n,Jm 


(6.15b) 


where  en  is  a  bound  state  eigenvalue  of  the  isotropic  component  v°  (z)  and 


s  + 


(6.16) 


en,Jm  =  1  <Jm|Pjl(cose)|Jm><njv^o(z)jn> 


(6.17) 


Thus  c^jm  are  the  zeroth  order  bound  vibration-rigid  rotor  (J3VRR)  eigenvalues 
and  ^ jm  are  the  first  order  perturbation  corrections  to  these  induced  by  the 
anisotropy  of  VQ0(z,  0). 

Pure  DSA  and  mixed  DSA-RMSA  resonances  are  obtained  from  Eq.  (.6.11)  with 
G  t  (00).  Making  the  diagonal  approximation  with  respect  to  J  again,  we  derive 
the  first  order  approximation  to  these  resonance  energies  as 

S.GJm  =  Hm  k?  +1!t  Ji^Ji  +  ^ 


(6.18) 


(6.9) 


y  (z,  0,  4>,  R)  =  E  Xejm(z)  Yjm(0» 

Gum 

the  following  set  of  close  coupled  equations  is  derived 


[77  +  ki  +  T  VJ1  ♦  »  -  <5,  +  5>2  -  T  J'J  +  »J 

az 


pu 

G.^.^2  <JmlVG-G,(z*  0’  J'm'  >  xG'JV(z) 


(6.10) 


with  the  initial  rotational  state,  and  I  the  moment  of  inertia  of  the  rigid 
rotor.  For  weakly  corrugated  surfaces  we  neglect  the  off  diagonal  G  coupling, 
retaining  only  the  rotational  coupling: 


+  1)  -  (K, 


+  e>2  - "  J(J  +  »wi) 


(6.11) 


■  J  <J"lV00tz-  6)lJ'",>XW'nl'lzl 

J  m  -Ti  ~ 

Since  the  initial  state  has  G  =  (.00)  we  may  drop  the  subscript  G  and  assume 

G  =  (00)  in  Eq.  (6.11).  The  possibility  of  pure  RMSA  thereby  arises,  i.e.  resonances 

observed  at  eigenvalues  of  Eq.  (6.11)  with  G  =  (00)  cause  effective  coupling  of 

the  incident  beam  (usually  =  0)  to  higher  J  states,  which  correspond  to  closed 

channels  asymptotically  at  the  value  of  kz  .  The  strong  rotational  inelasticity 

seen  for  the  asymmetric  rotor  HD  makes  this  an  ideal  case  in  which  to  see  pure  RMSA 

[50],  and  it  has  been  previously  shown  that  for  HD  on  smooth  metal  surfaces  (Pt(lll) 

[27]  and  Ag (111)  [28])  these  resonance  states  are  very  well  approximated  by  the 

first  order  approximation  to  the  eigenvalues  of  Eq.  (6.11)  with  G  =  (00).  We  shall 

denote  the  laterally  averaged  interaction  potential  for  the  heteronuclear  species  HD  by 

V00(z»  6)  in  order  to  distinguish  it  from  that  for  the  homonuclear  species  H2 

and  D?,  denoted  by  V  (z,  0).  V  (z,  6)  is  expanded  as 
c.  00  00 


approximation  of  free  atomic  motion  parallel  to  the  surface  in  the  resonant 
state 

-K2k2/2M  =  4i2  (K.  +  G)2/2M  +  en  (6.7) 


when  ep  is  a  bound  state  eigenvalue  of  the  isotropic  laterally  averaged 
potential  V0Q(z).  Deviations  from  this  free  particle  approximation  will  occur 
at  degeneracies,  when  there  is  strong  coupling  by  one  of  the  periodic  compon¬ 
ents  VG(z).  It  is  then  necessary  to  calculate  the  two  dimensional  band  struc¬ 
ture  of  the  adsorbed  particle  in  order  to  predict  the  positions  of  interfering  , 
resonances.  This  has  been  used  to  extract  information  on  t/>e  magnitude  of  the 
periodic  potential  components  Vg(z)  [  24,  48l. 

Equations  (6.1)  and  (6.7)  have  been  applied  to  the  scattering  of  H2  and 
D2  from  dielectrics  and  graphite  [49,  18].  Due  to  the  spherical  nature  of  the 
homonuclear  hydrogen  isotopes  and  the  low  rotational  excitation  probabilities  it 
is  to  be  expected  that  condition  (6.7),  which  neglects  the  molecular  degrees  of 
freedom,  will  also  be  valid  for  H 2  and  D2  to  a  first  approximation.  In  fact 
the  fhee  atom  approximation  will  be  even  better  for  scattering  from  metals  than 
from  dielectrics.  But  in  general,  the  internal  molecular  degrees  of  freedom 
will  affect  the  resonance  conditions.  This  will  now  be  illustrated  for  a  rigid 
rotating  diatomic  incident  on  a  corrugated  surface. 

The  generalization  of  Eq.  (1.1)  to  a  diatomic  molecule-periodic  surface 


system  is 


V(z,  9,  4>,  R)  =  E  V_( z,  9,  <j>)  e 
G  ~ 


(6.8) 


where  the  coordinate  system  is  the  standard  laboratory  fixed  frame,  and  a  rigid 
rotor  has  been  assumed  so  that  the  dependence  on  r  is  dropped.  For  a  flat  sur¬ 


face  the  laterally  averaged  term  must  be  independent  of  <{>,  i  .e.  VQ0(z,  6). 

When  the  molecular  wave  function  is  expanded  in  a  direct  product  basis  of  rota¬ 


tional  and  diffraction  states 


Born-Oppenheimer  formalism  proposed  by  Levi  and  Suhl  for  slow  collisions  with 
soft  crystals  [40]  and  deserves  further  attention. 


VI.  Selective  Adsorption  Kinematics:  The  Molecular  Case 

The  kinematics  of  selective  adsorption  are  determined  by  two  conditions; 
i)  conservation  of  parallel  momentum 


Kf  -  ^  +  6  (6.1) 

and  ii)  conservation  of  total  energy 

Ef  =  E1nt  (6-2) 

where  *  (K,  kz)  is  the  incident  particle  wavevector  with  components  K  =  (kx» 
ky)  parallel  and  kz  perpendicular  to  the  surface,  respectively,  M  is  the  total 
mass  of  the  incident  particle,  E^nt  its  internal  energy,  and  6  a  reciprocal  lat¬ 
tice  vector  of  the  surface. 

For  atom-solid  scattering,  with  the  interaction  potential  expanded  as  the 
Fourier  sum,Eq.  (.1.1), the  resonance  energies  are  determined  by  eigenvalues  of 
the  close  coupled  equations: 

Tt  *  A  -  <?i+?>2  XsW-  :  <t^l  (6.3) 


where 


Vz> 


(6.4) 


and  the  atomic  wave  function  has  been  expanded  in  a  plane  wave  basis 

*(r)  -  Z  XG(z)$G  (R) 

G  ~ 


(6.5) 


i(K.  +  G)  *R 
$G  W  =  As  e  *  ” 


(6.6) 


and  A„  the  area  of  the  surface  unit  cell.  In  most  cases  it  is  sufficient  to  retai 
s 

only  the  diagonal  term  on  the  right  hand  side  of  Eq.  (6.4),  giving  the  usual 
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fitted  slopes  of  linear  log  (I/IQ)  vs.  T&  plots  for  a  range  of  0  -*■  1  and 
0  -*•  2  rotationally  inelastic  peaks  measured  at  values  of  Ez  not  equal  to 
RMSA  resonances.  These  data  now  do  give  a  consistent  value  of  well  depth  D 
and  Debye  temperature  6p,  also  listed  in  Table  3.  However,  the  well  depth  of 
31.9  meV  is  considerably  smaller  than  the  value  of  45  meV  extracted  from  the  homo- 
nuclear  isotopes.  In  Figure  9  we  show  the  theoretical  Debye-Waller  correction  fac¬ 
tors  I/IQ  for  each  rotationally  inelastic  transition  as  a  function  of  Ez>  the  inci¬ 
dent  energy  normal  to  the  surface,  calculated  assuming  a  well  depth  of  a)  45.2  meV  and 
b)  31.9  meV.  Superimposed  on  these  lines  are  the  experimentally  measured  Debye- 
Waller  correction  factors  for  all  transitions  J  =  0  -*•  J  =  0,  1,  2.  We  note  that 
the  experimental  points  fall  within  the  envelope  of  the  calculated  values  for 
b),  but  not  for  a).  Thus  although  the  0  ■+■  0  and  0  -►  3  peaks  have  an  irregular 
behavior  with  energy  Ez  ,  they  do  always  fall  within  the  range  allowed  by  the 
total  rotational  inelasticity  at  any  given  Ez-  The  disagreement  of  the  curves 
with  the  experimental  values  in  the  upper  panel  of  Fig.  9  is  strong  evidence 
against  a  well  depth  of  45  meV.  Nevertheless,  the  discrepancy  between  the  internally 
consistent  and  D^  results  of  D  =  45  meV  and  the  HD  result  of  31  meV  leaves 
open  an  interesting  and  important  question  for  Debye-Waller  analysis.  We  also 
note  that  a  similar  type  of  result  was  seen  for  molecular  hydrogen  scattering 
from  Pt(lll)  [47],  where  the  H2  0  -*■  0  data  gave  a  well  depth  of  72  meV  while  the 
HD  rotationally  inelastic  data  gave  a  well  depth  of  55  meV. 

Our  conclusion  is  that  simple  Debye-Waller  theory,  Eq.  (5.1),  appears  to 
work  well  for  H^  and  D^  and  also  for  HD,  if  proper  allowance  is  made  for  the 
rotational  inelasticity,  but  that  the  results  from  the  homonuclear  and  hetero- 
nuclear  isotopes  are  quite  different.  The  anisotropic  component  of  the  laterally 
averaged  molecule-surface  potential  is  apparently  playing  an  important  role 
here,  which  has  yet  to  be  determined.  This  is  suggested  by  the  adiabatic 
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diffuse  background  scattering  has  been  subtracted  from  the  specular  peak.  The 

plots  show  very  good  linearity  over  a  range  of  temperatures,  showing  essentially 

only  the  expected  deviations  at  high  temperatures  due  to  anharmonic  contributions 

to  lattice  vibrations  (low  incident  angle  plots)  and  at  low  temperatures  the 

deviations  characteristic  of  zero  point  vibrational  energy  and  contamination 

effects  [42]  (high  incident  angle  plots).  Extrapolation  of  the  linear  segments 

to  0  K  gives  close  to  unit  reflectivities  (.83%  for  H^,  88%  for  D2,  94%  for  He) 

indicating  a  high  degree  of  surface  perfection.  Figure  6  presents  the  slopes 

for  each  of  the  least  squares  fitted  lines  of  Figure  5,  plotted  as  a  function 

of  E2  of  the  incident  beam.  Equation  (5.2)  is  then  used  to  obtain  an  estimate 

of  the  well  depth,  D,  and  surface  Debye  temperature  6p.  Calculated  mean  square 
2 

displacements  <uz>  at  273  K  are  also  listed.  These  values  are  reported  in 
Table  3  where  it  is  apparent  that  the  homonuclear  isotopes  H2  and  D2  give  a 
consistent  result  of  D  =  45  meV,  but  slightly  different  surface  Debye  temperatures. 
The  latter  are  greater  than  or  equal  to  the  bulk  value,  contrary  to  expectations 
based  on  the  relative  amplitudes  of  surface  and  bulk  atom  motions. 

Figure  7  shows  the  temperature  dependence  of  some  of  the  HD  rotational ly 
elastic  and  rotationally  inelastic  specular  (all  G  =  (00))  peaks  after  subtraction 
of  the  diffuse  background.  Efforts  were  made  touseEq.  (5.3)  to  extract  D  and  0D 
from  the  temperature  dependence  of  the  0  +  0  peak  in  the  angular  and  energy 
ranges  where  all  rotationally  inelastic  channels  are  closed  (E.  <  11  meV) .  No 
consistent  results  could  be  obtained,  either  coincident  with,or  between  RMSA 
resonances.  This  is  not  surprising  since  the  resonances  in  this  region  are  very 
strong  and  closely  spaced  in  energy.  Equation  (5.3)  is  only 
applicable  in  the  region  between  resonances,  which  is  ill-defined  for  the  0  +  0 
scattering  (Fig.  7).  When  rotationally  inelastic  channels  are  open  the  0  +  0 
peak  intensities  are  much  smaller,  and  are  subject  to  large  experimental  uncer¬ 
tainties,  as  are  the  0  +  3  intensities.  Figure  8  then  shows  the  least  squares 


where  and  6..  are  related  for  diffractive  scattering  by 


Eicos26f  =  E.  -^2(Ki  +  G)2/2M 


(5.4) 


and  by 


E.cos20f  =  E.cos2e.  +^2[J.(J.  +  1)  -  Jf(Jf  +  1)]/2I  (5.5) 


for  rotational ly  inelastic  scattering.  M  and  I  are  the  mass  and  moment  of  inertia  of 
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the  gas  particle,  M$  the  crystal  atom  mass,  D  the  well  depth  and  0p  the  surface  Debye 
temperature.  The  temperature  dependence  of  Eq.  (5.2)  is  often  used  to  extract  an  esti¬ 
mate  of  the  well  depth,  D.  Criticisms  of  this  form  include  the  ambiguity  in  the  defi¬ 
nition  of  Ak,  i.e.  whether  this  should  be  measured  at  the  potential  minimum  or  in  the 
repulsive  region  [39],  and  the  assumption  that  the  incoming  particle  collides 
with  only  one  surface  atom,  [41,  44]  which  is  certainly  not  the  case  for  large 
species  such  as  Ne  and  H2  [41,  45].  The  question  of  the  validity  of  Eq.  (5.1) 
itself  for  atomic  scattering  was  investigated  by  Levi  and  Suhl  [40],  who  showed 
that  for  fast  collisions  the  equation  does  hold  under  either  completely  coherent 
or  completely  incoherent  lattice  motion  and  the  assumption  of  pairwise  interactions. 
Thus  Debye-Waller  theory  is  applicable  to  both  neutrons  and  electrons  by  virtue 
of  the  short  collision  times  involved.  However,  for  the  majority  of  atomic  and 
molecular  scattering,  the  collision  time  is  of  the  same  order  as  the  characteristic 
vibration  time  of  the  crystal.  According  to  the  criterion  of  Levi  and  Suhl  [40], 
the  scattering  of  H2  from  Ag(lll)  is  actually  a  borderline  case. 

With  the  above  caveat  concerning  the  Debye-Waller  factor  in  mind,  but  lack¬ 
ing  a  more  suitable  alternative,  we  have  analyzed  the  temperature  dependence  of 
the  specular  rotationally  elastic  scattering  of  H2,  D2  and  rotationally  inelastic 
scattering  of  HD  from  Ag(lll)  with  the  simplest  model,  Eqs .  (5. 2-5. 3).  Figure  5 
shows  the  linear  dependence  of  the  logarithms  of  the  specular  reflectivities 
versus  surface  temperature  at  various  incident  angles  for  H0,  D2,  and  He.  The 


used  in  X-ray  scattering  [37],  and  more  recently  in  neutron  scattering  and 
in  LEED  [38] 


I/I 


o 


% 


2 

exp(-<(u*Ak)  > 

2  2 
exp(-  <uz>Ak  ) 


(5.1) 


Here  I  is  the  experimentally  measured  intensity  at  a  fixed  surface  temperature 

Ts,  IQ  the  hypothetical  intensity  for  a  stationary  lattice  at  0  k,  Ak  =  k^-k^ 

2 

the  wavevector  change  of  the  atom  upon  collision  with  the  surface  and  <u  > 
the  mean  square  displacement  of  a  surface  atom  at  Ts.  There  has  been  much 
controversy  over  the  applicability  of  such  a  Debye-Waller  factor  to  atomic 
and  molecular  scattering  from  surfaces  [39-45].  The  simplest  approach  is 
to  evaluate  Ak  at  the  potential  minimum,  thereby  including  the  effect  of  a 
thermally  averaged  attractive  well  (Beeby  correction  [46])  and  to  evaluate 
<u  >  assuming  i)  collision  with  a  single  surface  atom  ii)  <u  >  is  dominated 
by  its  normal  component  iii)  harmonic  lattice  motion  and  tv)  a  Debye  model 
for  the  solid  so  that  <uS  =  3-fi  Ts/MskB0Q-  The  resulting  expression  for 
specular  scattering  is 

I/IQ  =  exp[-  24  Mg(E.cos2e.  +  D)  yM^e2]  (5.2) 

and  for  diffractive  or  rotationally  inelastic  scattering 
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I 

I 

I 

[  were  ca.  8%  for  He  and  5%  for  using  300  K  nozzle  expansions.  They  correspond 

to  Debye-Waller  corrected  reflectivities  (i.e.  extrapolated  to  Ts  =  0  K)  of 
greater  than  80%. 

The  azimuthal  orientation  of  the  crystal  was  calibrated  by  optimizing  hydro¬ 
gen  diffraction  along  the  <112>  and  <101>  symmetry  directions.  Azimuthal  set¬ 
tings  falling  between  these  two  angles  (0°  and  30°)  were  set  using  a  pointer 
and  scale  attached  directly  to  the  crystal  mount.  This  pointer  could  be  set  to 
better  than  0.2°  by  using  a  remote  telescope  to  read  the  angle  scale.  Polar 
detector  and  crystal  angles  were  set  to  better  than  0.05°. 

In  each  experimental  session,  beam  characterizing  time-of-fl ight  measure¬ 
ments  were  taken  before  and  after  any  beam  scattering  measurements.  Three  dif¬ 
ferent  beam  energies  were  used  for  each  isotope,  shown  in  Table  2.  p-hydrogen 
and  o-deuterium  were  only  available  for  the  beams  of  lowest  energies  (nozzle 
size  12y) .  Two  types  of  scattering  data  were  collected:  1)  angular  distributions 
obtained  by  scanning  the  detector  at  0.2°  increments  for  fixed  incident  angle  and 
2)  peak  intensity  scans  obtained  by  measuring  the  maximum  intensity  of  a  given 
diffraction  peak  after  incrementing  the  incident  angle  by  ca.  0.2°.  In  practice, 
DSA  data  was  collected  by  alternately  moving  the  crystal  polar  angle  by  ca.  0.2°, 
and  then  scanning  the  detector  through  the  specular  peak  with  the  computer 
recording  in  real  time  the  peak  intensity  and  corresponding  incident  and  scatter¬ 
ing  angles  with  respect  to  the  surface  normal.  The  peak  heights  determined  in 
this  manner  include  contributions  from  both  coherently  scattered  particles  and 
the  diffuse  inelastic  background. 

V.  Debye-Waller  Analysis 

Inelastic  attenuation  of  elastic  atom-surface  scattering  intensities  by 
thermal  vibrations  and  zero  point  motion  of  the  surface  is  usually  analyzed  by 
rescaling  with  a  Debye-Waller  attenuation  factor  analogous  to  that  originally 
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0-D2  to  that  of  n-D^.  Note  also  the  disappearance  of  the  J  =  1  -+•  3  peak  for 
0^  after  conversion.  Thus  the  successful  generation  of  a  pure  J  =  0  beam  is 
confirmed.  The  generation  of  beams  richer  in  ground  rotational  state  molecules 
can  also  be  directly  demonstrated  by  detecting  a  small  increase  in  beam  transla¬ 
tional  energy  after  conversion  due  to  the  extra  rotational  energy  relaxing  to 
the  ‘  inslational  degree  of  freedom,  as  shown  for  hydrogen  in  Figure  4. 

IV.  Experimental  Procedure 

The  single  crystal  surface  Ag(lll)  was  spark-cut  from  a  99.999%  Ag  boule 
and  initially  polished  at  the  Cornell  Materials  Preparation  Laboratory.  It  was 
subsequently  repolished  in  our  laboratory  with  very  fine  cerium  oxide  powder, 
and  electro-chemically  etched  with  a  saturated  chromic  acid- 10%  HC1  solution  [36] 
to  remove  the  damaged  topmost  layers  which  result  from  mechanical  polishing.  The 
oxide  layer  formed  from  this  chemical  treatment  was  removed  with  a  10%  HNO^  solu¬ 
tion.  The  orientation  of  the  crystal  was  checked  by  Laue  X-ray  back-reflection, 
and  was  found  to  be  within  0.5°  of  the  (ill)  plane.  The  underside  of  the  crystal 
disc  was  spotwelded  to  two  1mm  thick  Pt  wires,  which  were  attached  to  the  manipu¬ 
lator  through  two  tungsten  rods.  A  Chromel-Alumel  thermocouple  was  also  spotwelded 
to  the  back  of  the  crystal.  The  crystal  was  resistively  heated  or  conductively 
cooled  through  the  Pt  support  wires.  Temperature  regulation  was  achieved  with  a 
home  built  circuit  which  controls  the  current  source. 

After  the  initial  bake-out,  the  crystal  was  subjected  to  several  cycles  of 
Ar+  sputtering  (2  uA  at  500  V  for  10  min.)  until  no  appreciable  contamination 
could  be  detected  with  a  single-pass  cylindrical  mirror  analyzer  Auger  spectro¬ 
meter.  The  crystal  was  then  annealed  at  600  K  for  one  hour  to  remove  sputter 
damage.  The  success  of  this  annealing  procedure  was  verified  by  the  fairly  uni¬ 
form  He  and  reflectivities  which  were  measured  at  several  points  along  the 
crystal.  Typical  reflectivities  at  6^  =  50°  with  the  surface  at  room  temperature 


The  values  of  the  experimental  bound  state  energies  for  p-H2»  n-H^.  o-D2, 
and  n-D2  are  listed  in  Table  4.  The  small  energy  differences,  AEn,  arising 
between  the  n-H^  (n-D^)  and  p-H2  (o-D2)  results  are  due  to  the  slight  shifts 
of  the  J  ^  0  m  f  0  resonance  energies  relative  to  |0,m>  =  |0,0>.  The  normal 
hydrogen  (deuterium)  results  represent  a  weighted  superposition  of 
the  |J,m>  =  |0,0>,  |l,0>,  and  |1,±1>  resonances.  Although  the  AEn's  are  just 
somewhat  larger  than  the  standard  deviations,  they  are  not  artifacts  due  to 
systematic  drifts.  This  was  confirmed  by  taking  data  for  p-H2  and  n-H^  sequen¬ 
tially,  without  changing  any  other  experimental  conditions.  This  was  accom¬ 
plished  utilizing  the  bypass  value  that  was  intentionally  built  into  the  p-H2 
(o-D^)  generator.  Note  that  since  n-D2  contains  only  33%  J  =  1  the  level 
shifts  between  o-D^  and  n-D^  are  smaller  than  those  between  p-H^  and  n-H^.  A 
more  complete  discussion  of  these  energy  shifts  is  presented  in  our  paper  dealing 
with  the  anisotropic  component  of  the  laterally  averaged  hydrogen/Ag( 111)  poten¬ 
tial  [13]. 


VIII.  Analysis  of  Bound  State  Data  From  DSA  Resonances:  p-H2/o-D2 
A.  Assignment  of  Quantum  Numbers 

For  o-D2  and  p-H2  the  J,m  dependent  shifts  of  the  bound  state  levels  are 
equal  to  zero  and  the  DSA  resonance  positions  provide  a  direct  measurement  of 
the  eigenvalues  of  the  isotropic  component  v°Q(z)  of  the  laterally  averaged  poten¬ 
tial  Eqs.  (6.16  and  6.18  with  J  =  J..  =  0).  At  the  present  time  there  are  two 
ways  to  generate  a  one  dimensional  gas-solid  potential  v°Q(z)  from  the  spectrum 
of  its  eigenvalues  {en>:  either  by  fitting  any  one  of  several  model  potentials, 
or  by  using  a  modified  RKR  procedure  to  generate  a  set  of  turning  points  together 
with  assumption  of  the  asymptotic  form  v°Q(z)^  -  C3/(z-zq)3,  Eq.  (1.2)  [30].  How¬ 
ever,  in  each  case  it  is  necessary  to  first  have  unambiguous  quantum  number  assign¬ 
ments  for  the  levels  {en)«  This  is  facilitated  by  independent  knowledge  of  the 


the  well  depth.  When  however,  as  is  the  case  here,  it  is  not  clear  that 
standard  Debye-Waller  analysis  is  valid  and  isotopic  Debye-Waller  studies  give 
inconsistent  results  for  the  well  depth  (.Section  V),  quantum  assignments  and  a 
crude  estimate  of  the  well  depth  may  usually  be  obtained  by  simultaneous  fitting 
of  different  isotopic  data  for  the  higher  vibrational  levels  to  the  asymptotic 
expression  of  Le  Roy  [30] 

lenl 1/6  =  2.02176  C"1/3  (nd-n)  (8.1) 

where  the  mass  reduced  vibrational  quantum  number  is  n  =  (v  +  and  is 

the  value  of  n  at  the  dissociation  limit. 

Figure  17  shows  the  Le  Roy  plot  for  four  different  sets  of  simultaneous 
assignments  for  the  and  levels.  We  note  that  with  the  exception  of  the 
lowest  bound  state  of  hydrogen,  the  plots  are  highly  linear  over  the  entire 
range  of  measured  bound  states  with  equal  slopes  for  both  species,  indicating 
that  Eq.  (8.1)  is  valid  over  the  measured  range  of  energies.  The  lowest  vibra¬ 
tional  level  for  H2  drops  below  the  limiting  slope.  In  each  case  best  straight 
lines  are  fit  separately  for  the  two  isotopic  species  to  test  their  closeness. 
The  gradients  and  intercepts  are  shown  in  Table  5.  Assignments  (a)  and  (d)  are 
both  very  close,  while  (b)  and  (c)  clearly  do  not  give  simultaneous  fits  to 
Eq.  (8.1)  for  both  isotopes,  and  so  can  be  rejected.  Any  other  possibilities 
show  poorer  correlation  than  these.  The  assignments  (a)  and  (d)  give  very  dif¬ 
ferent  positions  for  the  well  depth  which  is  well  illustrated  by  the  behavior 
for  the  bound  state  of  both  isotopes  at  'v  -10  meV.  For  (a)  this  H2  level  has 
n  =  2  and  D2  level  n  =  3.  This  assignment  gives  a  well  depth  of  %  32  meV  from 
the  model  potential  fittings  shown  later.  The  Y-intercept  of  the  fitted  lines 
in  (a)  provides  an  upper  bound  of  'v,  37  meV  for  the  well  depth,  and  implies  one 
unobserved  deeper  level  for  D?  and  none  for  H?.  The  slope  gives  a  crude  estimate 


of  C^j  which  we  obtain  as  983  ±  37  meV-A^  from  (a).  For  (d)  the  Hg  level 
at  'v  -10  meV  has  n  =  4  and  the  D2  level  n  =  6.  This  implies  two  unobserved 
deeper  levels  for  and  four  for  82*  and  yields  a  well  depth  of  ^  60  meV. 

It  is  clear  that  this  is  a  rather  unusual  situation  in  that  the  simul¬ 
taneous  isotopic  fitting  procedure  is  not  on  its  own  sufficient  to  disting¬ 
uish  between  two  very  different  quantum  number  assignments.  This  is  a  con- 

L2 )  where  n, 
u2 

and  n2  are  both  small  integers  Cn j  =  3,  n2  =  2).  Fortunately,  we  can  now 
use  the  rotationally  inelastic  transition  probabilities  for  HD  to  give  a  crude 
but  independent  estimate  of  the  well  depth,  which  is  sufficient  to  distinguish 
between  the  values  of  30  and  60  meV.  Scattering  calculations  show  that  the 
gross  features,  in  particular  the  energy  dependence  of  the  relative  rotational 
transition  probabilities,  are  not  very  sensitive  to  the  anisotropy  of  the  po¬ 
tential  V^z,  0)  for  HD,  but  are  strongly  affected  by  large  variations  in 
the  well  depth  [13,  47  ].  In  the  next  section  we  present  several  model  poten¬ 
tials  used  in  the  calculations  and  then  compare  the  results  of  these  quantum 
scattering  calculations  using  potential  parameters  obtained  by  least  square 
fits  to  assignments  (a)  and  (d),  to  the  experimental  rotationally  inelastic 
transition  probabilities  for  HD. 

B.  The  Isotropic  Potential  v°  (z) 

_ 00 

The  following  four  model  potentials  were  used  to  fit  the  bound  states  en 
of  v°  (z)  for  H~  and  D9: 
i)  Morse  potential 

v°QU)  =  D  {exp[-2cx(z-ze)]  -2exp[-a(z-ze)]  }  (8.2) 

where  D  is  the  well  depth,  a  the  inverse  range  parameter  and  zg  the  equili- 


sequence  of  the  mass  ratios,  in  particular  that  nj()^2)  ^  n ^ ( 


ii)  Lennard- Jones  9-3  potential 


voo(z)  =  33/2  0/2  {  C°/<z+2o>:)9  -Co/Cz+z0))]3  }  (.8.3) 

where  D  is  the  well  depth  and  X  the  range  parameter.  The  potential  minimum 
is  given  by  zfi  =  3^cr  -  zQ. 

iii)  Variable  exponent  potential  (VEP)  [52] 

v°Q(z)  -  D{  [1  +  X/ptz-ze)r2p  -  2  [1  +  X/p(z-ze)]"p}  (8.4) 

where  D  is  the  well  depth,  X  is  a  inverse  range  parameter,  -1  <  1/p  <  1,  and 
zg  is  the  equilibrium  position. 

iv)  Exponential-3  potential  [53] 

v°0(z)  =  A  e"bz  -  f ( z ) C/ z 3 
bz 

A  =  3  De  e/(bze-3) 

4  6  (8.5) 

C  =  bz*  D/(bze-3) 

-2z/a 

f ( z)  =  l-(2z/a0(l  +  z/ao)  +  1)  e  0 


where  D  and  zg  are  the  well  depth  and  equilibrium  position  when  f(z)  =  1, 
and  b  is  a  range  parameter.  For  f(z)  4  1  these  values  and  eigenvalues  of 
the  potential  are  determined  numerically  using  a  Numerov-Cooley  integrator. 
f(z)  is  a  cut-off  function  for  the  z  term  designed  to  go  smoothly  from 
f(z)  =  1  at  some  large  z  to  f (z )  =  0  at  z  =  0  [54]. 

The  Morse  potential  has  a  simple  form  and  analytical  eigenvalues  but 

it  decays  too  fast  in  the  long  range  attractive  region.  The  9-3  potential 

-3  -9 

shows  the  correct  z  dependence  in  the  attractive  component  but  the  z 

repulsive  term  is  arbitrary.  Analytic  eigenvalues  have  been  obtained  in 

the  WKB  approximation  for  this  [55]  (accurate  to  1%).  Both  of  these  are 

2-parameter  potentials.  The  variable  exponent  potential  is  a  3-parameter 


generalization  of  the  Morse  potential  [52]  which  has  had  considerable 
success  in  fitting  atom-surface  and  molecule-surface  vibrational  spectra. 
Approximate  analytical  eigenvalues  have  also  been  derived  for  this  poten¬ 
tial  (also  good  to  1%)  [52].  Finally,  the  exponential-3  potential  is  a 
3-parameter  potential  which  displays  the  correct  asymptotic  form  expected 
from  theoretical  analysis  at  both  large  and  small  z.  This  is  the  most  real¬ 
istic  model  potential,  but  has  the  practical  disadvantage  that  no  known 
analytical  expression  exists  for  its  eigenvalues,  which  must  be  solved  for 
numerically. 

The  potential  parameters  are  determined  by  simultaneous  fitting  of  the 

p-H2  and  o-D2  bound  state  energies,  which  are  listed  in  Table  6.  Table  7 

displays  the  experimental  and  calculated  eigenvalues  for  the  four  fitted 

2 

potential  forms.  The  quantity  x  is  used  to  compare  the  quality  of  the  fit: 

x2  -  I  [(Ev-E^)/ov]2  (8.6) 

2 

with  values  of  x  calculated  separately  for  each  isotope  and  then  added  to 
give  an  overall  value  for  the  simultaneous  fit.  As  expected,  the  two  para¬ 
meter  potentials  give  poorer  fits,  with  the  Morse  potential  worst  of  all  and 
9-3  somewhat  better.  The  variable  exponent  potential  is  better  than  the  Morse 
potential  by  an  order  of  magnitude, and  gives  a  fit  comparable  to  that  of  the 
exponential-3  potential.  The  four  best  fit  potentials  are  shown  in  Figure  18. 

C.  Rotationally  Inelastic  Quantum  Scattering  Calculations  for  HD:  Effect 
of  Well  Depth 

Quantum  scattering  calculations  were  performed  for  interaction  of  HD 


with  a  flat  surface  using  the  laterally  averaged  interaction  potential  obtained 
by  transformation  of  the  symmetric  rotor-surface  potential,  Eq.  (6.21),  to 
the  offset  center  of  mass 
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where  z',  z  are  equal  to  the  normal  distance  from  the  surface  to  the  geometric 
center  and  center  of  mass  of  HD,  respectively,  and  6  is  the  offset  of  the 
center  of  mass  from  the  geometric  center.  This  results  in 

VQ0(z,  e)  =  v°^(z  +  6COS0)  +  v^0(z  +  6 cose)  P2(.cose)  (8.8) 

»  V®  (z)  +  Z  Vt  P  (cose)  (8.9) 

00  1=  0,1,...  00  z 

In  our  previous  calculations  for  HD  [27,  28]  the  simplest  possible  form 

2  2  o 

of  the  homonuclear  anisotropic  term  v^Cz)  was  used,i.e.  v^Q( z)  =  g  v°Q(z). 

Although  this  is  insufficient  for  a  complete  and  consistent  description  of  the 

rotational  cross  sections  for  all  these  isotopes,  H2,  D2,  and  HD,  this  crude 

form  is  adequate  to  investigate  the  difference  between  a  well  depth  of  30  meV 

and  one  of  60  meV  since  the  anisotropy  is  small. 

The  scattering  calculations  were  carried  out  using  the  Wigner  R-matrix, 

2 

with  an  L  basis  of  6  rotational  x  60  translational  functions  [56]. 

Figure  19  shows  the  rotational  transition  probabilities  for  assignments 
(a)  (well  depth  =  31.535  meV)  and  (.d)  (well  depth  =  59.941  meV)  calculated 
using  the  VEP  potential.  Also  included  are  the  deconvoluted  pointwise  experi¬ 
mental  transition  probabilities,  including  subtraction  of  the  diffuse  back¬ 
ground,  instrumental  and  energy  deconvolution,  and  Debye-Waller  corrections 
[13].  Since  the  latter  was  performed  at  incident  angles  incremented  by  5°, 
resonance  structures  were  not  resolved.  Experimental  values  have  been  scaled 
to  give  a  total  cross  section  of  unity.  Gross  features  of  the  calculated 
rotational  transition  probabilities  are  similar  for  all  four  model  potentials 
with  a  given  well  depth.  Figure  19  and  calculations  performed  at  intermediate 
well  depths  (e.g.  D  =  43  meV  [28])  show  that  two  features  sensitive  to  the 
well  depth  are  i)  the  energy  of  cross-over  for  the  *  0  -*■  J  =  1  and 

=  0  •+  J  =  2  transition  probabilities  and  ii)  the  magnitude  of  the  J..  =  0  - 
J  =  3  transition  probability.  It  is  clear  that  there  is  significantly  better 


agreement  with  the  experimental  results  when  D  =  31.535  meV  (Fig.  22b)  than  when 
D  =  59.941  meV  (Fig.  22a).  Thus  we  conclude  that  the  well  depth  is  ^30  meV  and 
that  assignment  (a)  is  the  correct  assignment  for  the  bound  levels  of  the  homo- 
nuclear  isotopes.  These  conclusions  are  fully  consistent  with  the  level  assign¬ 
ments  and  isotropic  (VEP)  potential  that  have  recently  been  determined  for  H 2/ 
Ag(110)  [57]. 


IX.  Corresponding  States  Analysis  for  Physical  Adsorption 

Vidali,  Cole  and  Klein  [31]  have  recently  proposed  a  principle  of  corres¬ 
ponding  states  for  the  laterally  averaged  physisorption  potential  which,  when 
taken  with  the  known  asymptotic  form  Eq.  (1.2)  gives  a  reduced  Bohr-Sommerfeld 


quantization  relation 

0(ej)  =  (n  +  1/2 ) b  (9.1) 

=  I  [enr-g(x)]%  dx  (9.2) 

b  =  *tt/[C31/3  D1/6(2M)>s]  (9.3) 

where  3  en/D  is  a  reduced  energy  and  g(z)  is  a  reduced  potential 

v£q(z)  =  Dg(z*)  (9.4) 

z*  =  (z-ze)  ( D/C3) 1/3  (9.5) 

This  has  been  used  by  Vidali  et  al .  with  a  reduced  exponential-3  potential 
g( z*)  =  (.3/(u-3))exp(-z*u/xm)  -  l/(z*  +  xm)3  (9.6) 

where  xm  =  (1-  3/u  f/3  to  fit  all  known  bound  state  resonance  data  for  hydrogen  and 
noble  gas  scattering  from  surfaces  to  plots  of  J(en)  vs.  en-  It  is  found  that 
all  alkali  halide  data  are  characterized  by  a  parameter  u  =  5.25  while  the  Hg/ 


noble  metals  and  noble  gases/graphite  systems  are  characterized  by  a  value  of 
u  =  6.5. 

The  success  of  this  principle  of  corresponding  states  reflects  the  domin¬ 
ance  of  the  attractive  component  of  v°Q(z)  in  determining  bound  state  energies. 
Since  from  the  high  linearity  of  the  LeRoy  plot.  Fig.  17,  this  is 


also  true  for  the  isotropic  component  v°Q(z)  of  the  laterally  averaged 
potential  for  molecular  hydrogen  on  Ag(lll),  the  bound  states  of  v°0(z)  are 
expected  to  comply  with  the  reduced  quantization  plots  obtained  by  Vi dal i 
et  al.  [31]  from  atomic  scattering  data.  Fig.  20.  shows  the  reduced  quanti¬ 
ties  J(en)  (Eq.  9-1)  plotted  against  en  for  each  of  the  quantum  number 
assignments  (a)  and  (d).  The  value  of  used  was  determined  from  the  slope 
of  the  Le  Roy  plots,  Eq.  (8.1).  The  solid  and  dashed  lines  represent  the 
reduced  plots  of  Vi  dal i  et  al .  obtained  using  Eqs .  (9.2)  and  (9.6)  with 
u  =  6.5  and  5.25,  respectively.  Assignment  (a)  is  in  good  agreement  with  the 
u  =  6.5  curve  appropriate  to  noble  metals,  while  points  from  assignment  (d) 
lie  somewhat  higher  than  the  u  =  5.25  curve  appropriate  to  alkali  halides, 
which  renders  this  assignment  less  plausible.  The  point  of  this  comparison 
is  not  to  derive  a  value  of  u,  which  is  specific  to  the  parametric  form, 

Eq.  (9.6),  but  to  match  the  bound  state  behavior  of  the  H2(D2)/Ag( 111)  system 
against  that  of  other  species  physisorbed  on  similar  surfaces  but  with  dif¬ 
ferent  well  depths.  Thus  it  is  of  special  interest  and  value  to  this  study 
because  of  the  particular  combination  of  isotopic  masses,  which  rendered  the 
results  of  an  isotopic  study  based  on  the  Le  Roy  analysis  inconclusive. 

Figure  21  shows  the  reduced  exponential-3  potential,  Eq.  (9.6),  for 
u  =  6.5  plotted  together  with  the  reduced  variable  exponent  potential,  Eq.  (8.4), 
derived  from  assignments  (a)  and  (d)  (see  Table  6).  The  attractive  portion 
of  the  VEP  potential  with  well  depth  D  =  31.535  meV  is  nearly  coincident  with 
that  of  the  reduced  exponential-3  potential,  as  expected  from  the  degree  of 
agreement  shown  in  Fig.  20. 

X.  Thermal  Desorption  Experiments 

Efforts  were  also  made  to  obtain  an  independent  estimate  of  the  hydrogen/ 
Ag(lll)  well  depth  from  thermal  desorption  and  modulated  beam  measurements. 
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The  experimental  pbservatiQn  was  that  no  substantial  adsorption  or  desorption 
of  Hg  (or  Dg)  was  detected  at  temperatures  >  19  K,  the  lowest  surface  tempera¬ 
ture  that  could  be  reached  with  our  crystal  mount.  Attempts  were  also  made  to 
examine  reflected  molecule  waveforms  from  a  19  K  surface;  no  noticeable  demo¬ 
dulation  was  observed.  These  facts  are  only  consistent  with  the  well  depth 
assignment  being  D  =  32  meV,  as  the  alternate  D  «  60  meV  well  depth  would  lead 
to  appreciable  hydrogen  adsorption.  This  can  be  demonstrated  using  the  first 
order  rate  equation  which,  for  a  constant  temperature  ramping  rate  B,  relates 
this  activation  energy  for  desorption  to  the  peak  desorption  temperature  [58]: 

Ea/RTp  =  (vl/8)  e><P(-Ea/RTp)  (10.1) 

where  Efl  is  the  activation  energy,  Tp  the  peak  desorption  temperature,  and 
v1  the  pre-exponential  factor  for  first  order  desorption.  This  equation  was 
solved  numerically  using  the  experimental  temperature  ramping  rate  of 
^1.3  K/sec  and  an  assumed  (typical)  frequency  factor,  =  10  sec  .  The 

calculated  peak  desorption  temperatures  which  result  from  this  exercise  are 
ca.  13  K  for  a  well  depth  of  32  meV  and  24  K  for  D  =  60  meV--  consistent  with 
our  earlier  conclusion  from  the  Le  Roy  analysis  and  rotationally  inelastic 
HD  scattering  that  the  lower  value  of  the  well  depth  is  the  correct  assign¬ 
ment. 


XI.  Summary  and  Conclusion 

In  this  paper  we  have  presented  a  comprehensive  study  of  the  spatially 
isotropic  component  of  the  laterally  averaged  molecular  hydrogen/Aq( 111) 
physisorption  potential.  Our  analysis  began  by  mapping  out  the  diffractive 
and  rotationally  mediated  selective  adsorption  scattering  resonances  for 
n-H2,  p-H2»  n-D2,  and  o-D2  in  the  Kx-Ky  reciprocal  lattice  plane.  This  was 
carried  out  as  a  function  of  incident  polar  angle  for  several  crystal  azimuths 
and  beam  energies.  The  resonances  show  small  energy  shifts  and  linewidth 
changes  between  p-H2  and  n-H2>  and  between  o-D2  and  n-D2,  respectively,  which 


are  a  direct  consequence  of  the  anisotropic  component  of  the  laterally  averaged 
molecule-surface  potential.  The  detected  differences  represent  an  average  over 
the  |J,m>  state  distributions  present  in  the  incident  n-h^  and  n-D^  beams.  The 
ability  to  resolve  DSA  dips  on  Ag(lll),  a  surface  of  extremely  low  corrugation, 
is  important  in  that  it  suggests  that  future  bound  state  studies  should  be  fea¬ 
sible  on  a  much  wider  selection  of  crystalline  surfaces  than  previously  thought 
possible.  Measurements  on  Pt(lll)  also  support  this  contention  [14]. 

The  experimentally  determined  bound  state  energies  for  o-D^  and  p-l^. 
i.e.  |J,m>  =  |0,0>,  have  been  used  to  obtain  the  unperturbed  bound  eigenvalues 
of  the  potential  well.  Best  fit  Lennard- Jones ,  Morse,  variable  exponent, and 
exDonential-3  potentials  having  well  depths  of  ^32  meV  are  derived  from  the 
data.  Unfortunately,  this  well  depth  assignment  could  not  be  made  exclusively 
on  the  basis  of  a  Le  Roy  type  analysis  [30],  as  two  possible  quantum  number 
assignments,  corresponding  to  well  depths  of  ca.  32  and  60  meV,  were  indicated. 
Comparison  of  experimental  HO  rotationally  inelastic  transition  probabilities 
with  results  from  theoretical  scattering  calculations  for  each  well  depth  in¬ 
dicates  that  the  true  well  depth  is  ^32  meV.  This  conclusion  is  further  sud- 
ported  by  thermal  desorption  measurements,  and  is  consistent  with  a  reduced 
("universal")  potential  which  has  been  formulated  for  closed-shel 1/noble  metal 
physisorption  systems  [31].  The  derived  isotropic  potential  is  also  in  excel¬ 
lent  agreement  with  recent  ab  initio  predictions  for  H2/Ag( 111)  [59]. 

Standard  Debye-Waller  analysis  fails  to  give  a  consistent  well  depth  for 
all  three  isotopes  H2,  O2  and  HD.  Since  this  system  involves  a  borderline 
fast/slow  collision  with  a  soft  lattice,  according  to  the  criteria  proposed 
by  Levi  and  Suhl  [40],  it  is  not  clear  that  the  standard  Debye-Waller  analysis 
should  be  valid.  The  significant  difference  between  the  heteronuclear  result 
and  that  for  the  homonuclear  isotopes  indicates  that  the  anisotropic  component 
of  the  potential  may  be  of  importance  to  thermal  attenuation  in  a  manner  which 
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is  not  well  understood  at  the  present  time. 

We  close  by  noting  that  the  isotropic  potential  determined  in  this 
paper  plays  an  essential  role  in  subsequent  studies  dealing  with  the 
spatially  anisotropic  [13]  and  periodic  [60]  component  of  the  hydroqen/Ag( 111) 
potential . 
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d2  converter 


p-H2  or  o-D2 


Figure  17 


Le  Roy  plots  of  the  p-^  (+)  and  o-D^  (•)  eigenvalues  for  four 
different  quantum  number  assignments:  a)  corresponds  to  D%32  meV; 
b)  one  extra  level  for  both  and  DgS  c)  one  extra  level  for  Hg 
and  two  extra  levels  for  D£;  d)  two  extra  levels  for  Hg  and 
three  extra  levels  for  D ^  corresponding  to  D  ^  60  meV. 

Figure  18  Plots  of  the  four  model  potentials  derived  from  the  experimentally 
determined  p-H^  and  o-D^  bound  state  energies  on  Aq (111) . 

Figure  19  Comparison  of  the  calculated  scattering  probabilities  for 

HD/Ag( 111)  based  on  D  =  59.9  meV  (upper)  and  D  =  31.5  meV  (middle) 
with  the  experimental  results  derived  from  angular  distributions 
collected  at  incident  angle  increments  of  5°. 

Figure  20  J(e)  defined  by  Eqs .(9.1-9.3)  versus  reduced  bound  state  energy 
for  p-H^  and  o-Dg/Agflll)  based  on  two  possible  well  depths. 

The  D  =  31.5  meV  values  fall  on  the  curve  appropriate  for  closed 
shell/noble  metal  adsorption  systems. 

Figure  21  Comparison  of  the  reduced  experimental  VEP  potentials  (based  on 

two  possible  well  depth  assignments)  with  the  reduced  corresponding 
states  exp-3  potential.  The  excellent  agreement  between  the 
D  =  31.5  meV  VEP  potential  and  the  corresponding  states  potential 
in  the  attractive  (relevant)  region  supports  the  contention  that 
this  is  the  correct  well  depth  for  hydrogen/Ag( 111) . 
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Figure  10 


Figure  11 


Figure  12 


Figure  13 


Figure  14 


Figure  15 


Figure  16 


Specular  peak  intensity  of  scattering  from  Aq( 111)  as  a 
function  of  incident  polar  angle  at  two  azimuths  symmetric 
about  the  <112>  direction.  The  presence  of  weak  but  reproducible 
dips  in  these  two  early  (low  signal-to-noise)  data  runs  confirmed 
that  diffractive  selective  adsorption  resonances  could  be  observed 
on  Ag(lll). 

Specular  peak  intensity  of  p-H2  scattering  from  Aq (111)  as  a 
function  of  incident  polar  angle  for  eight  successive  azimuthal 
angles . 

Selective  adsorption  loci  for  p-H2  plotted  in  the  surface  reci¬ 
procal-lattice  plane.  Arc  labels  denote  bound  levels  (meV)  and 
coupling  G  vectors.  Points  represent  dips  from  the  resonance 
spectra.  The  dashed  arc  denotes  the  energy  limit  for  the  p-H2 
beam. 

Specular  peak  intensity  of  o-D2  scattering  from  Ag(lll)  as  a 
function  of  incident  polar  angle  for  seven  successive  azimuthal 
angles . 

Selective  adsorption  loci  for  o-D2  plotted  in  the  surface  reci¬ 
procal-lattice  plane.  Arc  labels  denote  bound  levels  (meV)  and 
coupling  G  vectors.  Points  represent  dips  from  the  resonance 
spectra.  The  dashed  arc  denotes  the  energy  limit  for  the  o-D2 
beam. 

Specular  peak  intensity  of  n-H2  scattering  from  Ag(lll)  as  a 
function  of  incident  polar  angle  for  three  beam  energies. 

Specular  peak  intensity  of  n-D2  scattering  from  Ag( 111)  as  a 
function  of  incident  polar  angle  for  three  beam  energies. 
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FIGURE  CAPTIONS 


Figure  1  In-line  p-H2  and  o-D2  converter. 

Figure  2  Angular  distributions  of  n-H2  (upper)  and  p-H2  (lower)  scattering 

from  Ag(lll).  The  approximate  four- fold  increase  in  the  J  =  0  2 

transition  probability  in  the  lower  panel  confirms  that  the  gas 
has  been  successfully  converted  to  p-H2  (J  =  0  after  expansion). 

Figure  3  Angular  distributions  of  n-D2  (upper)  and  o-D2  (lower)  scattering 

from  Ag(lll).  The  approximate  50%  increase  in  the  J  =  0  -*■  2 
transition  probability  in  the  lower  panel  confirms  that  the  qas 
has  been  successfully  converted  to  o-D2  (J  =  0  after  expansion). 
Also  note  that  the  J  =  1  -►  3  transition  is  absent  in  the  lower 
panel . 

Figure  4  Time-of-flight  distributions  for  n-H2  and  p-H2*  The  higher  term¬ 

inal  velocity  of  the  p-H2  beam  confirms  that  the  p-H2  converter 
functions  properly. 

Figure  5  Debye-Waller  plots  for  H2>  D2,  and  He  specular  scattering  from 

Ag(lll). 

Figure  6  Linear  fits  to  the  slopes  of  the  lines  shown  in  Fig.  5  plotted 

versus  Ez- 

Figure  7  Debye-Waller  plots  for  various  HD  elastic  and  rotationally  in¬ 

elastic  peak  intensities  versus  T  . 

Figure  8  Linear  fits  to  the  slopes  of  the  lines  shown  in  Fig.  7  plotted 

versus  Ez- 

Figure  9  Theoretical  reciprocal  Debye-Waller  attenuation  plots  for  HD/ 

Ag(lll)  for  two  different  well  depths.  Superimposed  on  the 
calculated  curves  are  the  experimentally  determined  values.  These 
are  designated  by  numerals  which  indicate  the  final  rotational 
quantum  state  of  the  scattering  transition. 
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Tajhle  8 

Reduced  quantities  for  the  "universal  physisorption  potential". 


(a) 

(d) 

D 

= 

31.54  meV 

D 

=  59.94  meV 

C3 

= 

983  meV-A3 

C3 

,  =  829  meV-A3 

X 

= 

3.148  A 

1 

X 

=  2.400 

o 

A 

b  =  .1810 

b  =  .1722 

n 

en 

J  U„) 

n 

en 

J(V 

5 

.0384 

.996 

7 

.0202 

1.292 

4 

.0828 

.815 

6 

.0435 

1.119 

H2 

3 

.1737 

.634 

5 

.0913 

.947 

2 

.3205 

.451 

4 

.1687 

.775 

1 

.5349 

.272 

3 

.2814 

.603 

0 

.8161 

.091 

2 

.4294 

.431 

b  =  .1281 

b  =  .1218 

n 

en 

0  (en) 

n 

en 

J<v 

7 

.0358 

.961 

10 

.0189 

1.279 

6 

.0758 

.833 

9 

.0399 

1.157 

5 

.1259 

.705 

8 

.0662 

1.035 

D2 

4 

.2083 

.576 

7 

.1096 

.914 

3 

.3199 

.448 

6 

.1683 

.792 

2 

.4588 

.320 

5 

.2414 

.670 

. ...  ^  , .. ....  •■<  ...  i.  «L-  wmw  i.i «.  «j  A  TO  '.W.V.  ■. 1 
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Table  7 

Experimental  (En)  and  calculated  (en)  eigenvalues  of  the  Hydrogen/AgClll) 
physisorption  potential  for  four  model  potentials. 


n 

En 

en:  Morse 

9-3 

VEP 

Exp- 3 

5 

-1.21  ±  .06 

-.32 

-1.46 

-.97 

-1.09 

4 

-2.61  ±  .13 

-2.11 

-2.93 

-2.58 

-2.58 

H2  3 

-5.47  ±  .10 

-5.48 

-5.50 

-5.44 

-5.37 

2 

-10.11  ±  .09 

-10.42 

-9.70 

-9.98 

-9.93 

1 

-16.87  ±  .07 

-16.94 

-16.30 

-16.62 

-16.63 

0 

-25.74  ±  .07 

-25.03 

-26.27 

-25.81 

-25.67 

2 

X 

297 

139 

31 

22 

7 

-1.13  ±  .05 

-.54 

-1.68 

-1.20 

-1.30 

6 

-2.39  ±  .10 

-1.87 

-2.75 

-2.37 

-2.39 

5 

-3.97  ±  .06 

-3.98 

-4.33 

-4.15 

-4.10 

°2  4 

-6.57  ±  .10 

-6.87 

-6.62 

-6.69 

-6.60 

3 

-10.09  ±  .06 

-10.56 

-9.83 

-10.13 

-10.07 

2 

-14.47  ±  .20 

-15.04 

-14.24 

-14.62 

-14.60 

1 

— 

-20.30 

-20.19 

-20.33 

-20.30 

0 

— 

-26.36 

-28.09 

-27.41 

-27.20 

2 

X 

244 

192 

14 

17 

TOTAL  x2 

541 

331 

45 

39 

Table  6 
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Fit  parameters  for  the  model  potentials.  Those  listed  under  (a) 
with  D  ^  32  meV  are  the  correct  potentials  for  Hydrogen/Ag(lll) . 


Morse 

9-3 

(a) 

VEP 


D  [meV]  s,zg  [A]  a,X,B  [1/A] 

29.678  .8717 

32.908  2.4651 

31.535  1.0646  4.2920 
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Table  4 

Experimentally  determined  bound  state  energies  of  p-H2,  n-H2»  o-D2»  and 
n-02  on  Ag(lll).  a  is  the  standard  deviation,  m  the  number  of  data  points 
and  n  the  assigned  quantum  number.  AEp  is  the  difference  between  the 
energies  of  bound  p-H2  (o-D2)  and  (quantum  state  averaged)  n-H2  (n-D2). 


CSJ 

1 

CL 

E.  =  18.57 
b 

meV 

n-H2  Eb  = 

18.09  meV 

n 

En 

o  (m) 

En 

a  (m) 

AEn 

5 

-1.21 

.06(11) 

-1.12 

.05(5) 

.09 

4 

-2.61 

.13(22) 

-2.52 

.07(14) 

.09 

3 

-5.47 

.10(45) 

-5.35 

.09(30) 

.12 

2 

-10.11 

.09(42) 

-9.91 

.09(31) 

.20 

1 

-16.87 

.07(28) 

-16.55 

:il(14) 

.32 

-25.74 

.07(2) 

-25.53 

.21(2) 

.21 

o-02 

Eb  * 21 

.68  meV 

n 

En 

a  (m) 

En 

a  (m) 

AEn 

7 

-1.13 

.05(11) 

-1.17 

•12(3) 

-.04 

6 

-2.39 

.10(26) 

-2.36 

.06(6) 

.03 

5 

-3.97 

.06(21) 

-3.94 

.07(9) 

.03 

4 

-6.57 

.10(32) 

-6.52 

.07(10) 

.05 

3 

-10.09 

.06(25) 

-10.06 

•07(7) 

.03 

2 

-14.47 

.20(8) 

-14.40 

.05(4) 

.07 

2 

-15.25* 

.04(5) 

•  •  • 

.  •  . 

•  •  • 

1 

-20.97* 

.04(13) 

-20.98* 

.02(3) 

•  •  • 

*RMSA  dips  due  to  J  =  0  2  transitions 
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Table  3 

Results  of  Debye-Waller  Analysis.  D  is  the  derived  well  depth,  6p 

2 

the  surface  Debye  Temperature,  and  <uz>  the  mean  square  displacement 
at  273  K. 


Eb  [meV]  D  [meV]  6D  [K]  <u2>  [A2] 


E.  ImeV] 

D 

Pn  [psig] 

D  [y] 
n  L^J 

Av/v 

77.1 

200 

20 

5.0% 

32.3 

60 

20 

5.7% 

18.1 

100 

12 

4.5% 

82.8 

280 

20 

7.0% 

36.7 

80 

20 

7.2% 

20.7 


90 


12 


6.6% 


37 


Table  1 

Calculated  incident  beam  rotational  state  populations  for  n-H2, 
p-H2»  n-Dg,  and  o-D2  and  the  relative  probablities  of  the  J  =  0  2 

scattering  transitions.  Tr  is  the  rotational  temperature  of  the 
terminal  expansion.  and  W$  are  the  scattering  probabilities  of 
the  J  =  0  -*•  2  and  specular  peaks,  respectively.  Tq  =  302  K. 


[meV] 

Tr  [K] 

m 

■ 

m 

n 

Wp/Ws 

n-H2 

78.2 

143 

21.8 

74.5 

3.2 

0.5 

.0059 

p-h2 

85.3 

139 

88.2 

0.2 

i 

11.6 

0.0 

.0283 

81.1 

118 

42.5 

31.4 

24.0 

2.0 

.0527 

o-D2 

81.8 

116 

63.1 

2.1 

34.4 

0.1 

.0837 
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